Coulomb correlations in the honeycomb lattice: role of translation symmetry 
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The effect of Coulomb correlations in the half-filled Hubbard model of the honeycomb lattice is 
studied within the dynamical cluster approximation (DCA) combined with exact diagonalization 
(ED) and continuous-time quantum Monte Carlo (QMC). The important difference between this 
approach and the previously employed cluster dynamical mean field theory (CDMFT) is that DCA 
preserves the translation symmetry of the system, while CDMFT violates this symmetry. As the 
Dirac cones of the honeycomb lattice are the consequence of perfect long-range order, DCA yields 
semi-metallic behavior at small onsite Coulomb interactions U, whereas CDMFT gives rise to a spu- 
rious excitation gap even for very small U. This basic difference between the two cluster approaches 
is found regardless of whether ED or QMC is used as the impurity solver. At larger values of U, 
the lack of translation symmetry becomes less important, so that the CDMFT reveals a Mott gap, 
in qualitative agreement with large-scale QMC calculations. In contrast, the semi-metallic phase 
obtained in DCA persists even at U values where CDMFT and large-scale QMC consistently show 
Mott insulating behavior. 



I. INTRODUCTION 

The possible existence of a spin-liquid phase on the 
honeycomb lattice has recently attracted considerable at- 
tention. Meng et alX investigated the Hubbard model 
for this system at half-filling, using large-scale quantum 
Monte Carlo (QMC) calculations for clusters containing 
up to 648 sites. Careful finite-size extrapolations indi- 
cated semi-metallic behavior for onsite Coulomb interac- 
tions in the range U < 3.5i (t is the nearest neighbor hop- 
ping) and an antiferromagnetic insulator for U > 4.3i. 
The intermediate range 3.5t < U < 4.3i then corre- 
sponds to a Mott phase without long-range order, the 
hallmark of a spin liquid. These findings were, however, 
disputed by Sorella et al£ who performed similar QMC 
calculations for even larger clusters including up to 2592 
sites. The new results showed a considerably reduced 
spin-liquid phase, confined at most to the narrow win- 
dow 3.8t < U < 3M. 

The effect of nonlocal Coulomb correlations on the 
honeycomb lattice was also studied within the cluster 
extension of dynamical mean field theory^ (CDMFT). 
Wu et alA used continuous-time quantum Monte 
Carlo£ (CTQMC), whereas Liebsch- employed a multi- 
orbital-multi-site extension 7,8 of finite-temperature exact 
diagonalization- (ED) as impurity solver. Despite the 
fact that in the ED CDMFT calculations it was possible 
to include only a relatively small bath (six bath levels 
per six-site unit cell), the cluster self-energy components 
were found to be in nearly quantitative agreement with 
the CTQMC CDMFT results. (For a detailed compari- 
son see Fig. 25 of Ref.Hol.) In particular, for U ~ 5t both 
schemes revealed a Mott phase, with an excitation gap 
A w 0.6i, in close agreement with the one found by Meng 
et a/i With decreasing U, the CTQMC results at tem- 
peratures T > O.Obt indicated the closing of the Mott gap 
near U = 3.8i^ while the ED results at lower tempera- 
ture T — 0.005t revealed a weak insulating contribution 
to the self-energy at the Dirac points at arbitrarily low 



U & For U < 3t the small gap associated with this self- 
energy was, however, difficult to resolve in the spectral 
distributions due to the temperature rounding of the gap 
edges. 

Analogous ED CDMFT calculations (also for six bath 
levels) were carried out by He and Luii at a considerably 
lower effective temperature (T = 10~ 5 t). The excitation 
gap in this case was found to extend to U — >• 0. On the 
basis of these results the authors concluded that the spin- 
liquid phase of the honeycomb lattice at half-filling exists 
from U — up to the onset of the antiferromagnetic 
phase near U — 4.5i. 

Closely related to these works are two calculations 
based on the variational cluster approximation-^ (VCA) 
by Yu et alJ£ and Seki and OhtaJ^. In both cases, ED 
was used as impurity solver, with six bath levels as in 
Refs.iand[U Whereas Yu et al. identified a spin-liquid 
phase in the range U ~ 3t — 4t and semi-metallic be- 
havior at smaller values of U, Seki and Ohta obtained a 
similar insulating contribution to the self-energy at the 
Dirac points as in Ref. M and concluded that the Mott 
gap persists down to arbitrarily small values of U. 

Most recently, Hassan and Senechal 15 performed 
ED calculations for the honeycomb lattice within 
VCA, CDMFT and the cluster dynamical impurity 
approximation^ (CDIA). They argued that a bath con- 
sisting only of six levels is insufficient and leads to the 
erroneous conclusion that the system is gapped for all 
nonzero values of the onsite Coulomb interaction U. 

In this context it is also important to recall the results 
of functional renormalization group (FRG) calculations 17 
for the honeycomb lattice which reveal a stable semi- 
metallic phase below about U « 3.8£. 

In view of these contradictory results it is evident that 
the possible existence and extent of the semi-metallic 
phase of the honeycomb lattice are difficult to determine 
within present non-local many-body techniques. In par- 
ticular, it is not clear which assumptions and approxi- 
mations give rise to certain consequences: the size of the 
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correlated cluster, the size and symmetry of the bath in 
ED, the accessible temperature range, the accuracy of 
spectral functions at very low energies, etc. Naturally, 
these uncertainties also affect the identification of the 
elusive spin-liquid phase. 

The purpose of this work is to shed light on some of 
these issues by comparing new results derived within the 
dynamical cluster approximation 15 (DCA) with previous 
ones obtained within CDMFT As impurity solver we 
use finite-temperature ED as well as CTQMC. The nearly 
quantitative agreement between the ED and CTQMC 
self-energies, within DCA as well as CDMFT, demon- 
strates that the intrinsic limitations of these impurity 
solvers are not the cause of the discrepancies between 
the various results cited above. 

Instead we show here that, in the special case of the 
honeycomb lattice, it is of crucial importance to preserve 
the translational invariance of the system. Obviously, 
any deviation from bulk symmetry opens a gap at the 
Dirac points. Recall, for instance, the single-particle gaps 
obtained for graphene ribbons. Thus, the semi-metallic 
and spin liquid phases can only be studied properly by 
using many-body methods that do not violate translation 
symmetry. This argument disqualifies CDMFT which is 
well-known to yield a self-energy that is not translation- 
ally invariant j 18 ' 19 The self-energy components in this 
scheme account for correlations within the unit cell, but 
not between cells. We therefore believe that all CDMFT 
calculations performed until now for the honeycomb lat- 
tice should exhibit, at low U and low T, an excitation 
gap which is an artifact caused by the lack of translation 
symmetry of the self-energy. Although this gap is related 
to the presence of the local Coulomb interaction, it is not 
a true Mott gap but merely the consequence of the in- 
trinsic limitation of the cluster approach. As a result, 
CDMFT and other schemes that do not preserve trans- 
lation invariance are not suitable for the identification of 
a spin-liquid phase on the honeycomb lattice. 

The comparison of the CDMFT self-energy with anal- 
ogous results derived within DCA, for ED as well as 
CTQMC, underlines this point. In DCA, the self-energy 
is by construction translationally invariant, so that the 
electronic structure at low U is semi-metallic, in agree- 
ment with the predictions based on large-scale QMC and 
FRG calculations^^ The spurious tail of the excita- 
tion gap at small U and low T that is seen in CDMFT is 
absent in DCA. 

As will be shown below, in the case of the honey- 
comb lattice, the DCA condition that ensures translation 
symmetry is too rigid for the description of correlations 
within the unit cell. As a result, the semi-metallic phase 
is still stable near U — 5t — 6t where CDMFT and large- 
scale QMC calculations already find Mott insulating be- 
havior. Thus, CDMFT and DCA may be viewed as com- 
plementary cluster schemes: DCA is preferable at low U 
since it maintains the long-range order that is crucial for 
the Dirac cones, whereas CDMFT yields a more realistic 
description of short-range correlations in the Mott phase 



when the absence of translation symmetry plays a minor 
role. 

We also note here that the gap tail obtained in 
CDMFT at small U is not related to the finite size and 
symmetry properties of the bath used in ED. On the 
contrary, in the special case of the honeycomb lattice, a 
rather small bath containing only six levels is sufficient 
for the description of short-range correlations within the 
six-site unit cell. The reason is that, because of the semi- 
metallic properties of the system, the projection of the 
bath Green's function on a finite-cluster is not affected 
by the usual low-energy disparities that arise in the case 
of correlated metals. 

The outline of this paper is as follows: In Section II we 
discuss the application of DCA and CDMFT to the hon- 
eycomb lattice and point out the key difference between 
the self-energies obtained within these schemes. Section 
III presents the main ingredients of the ED impurity 
solver for both DCA and CDMFT. Section IV provides 
the discussion of the results obtained within ED DCA, 
and the comparison with analogous CTQMC DCA re- 
sults. The summary is presented in section V. Through- 
out this work only paramagnetic phases are discussed. 



II. DCA VS. CDMFT FOR THE HONEYCOMB 
LATTICE 

To describe Coulomb correlations in the honeycomb 
lattice we consider the single-band Hubbard Hamiltonian 

H = -t (4a c i° + H - c -) + U J2 n ^4, (!) 
(ij)a i 

where t is the nearest neighbor hopping term and U the 
on-site Coulomb energy. Throughout this paper t = 1 
defines the energy scale. The non-interacting band dis- 
persion is given by: e(k) = ±t\l+e lk:c ^+e l ( k:c ^ +ky3 ^ 2 \. 
The nearest neighbor spacing is assumed to be a = 1. 
We choose a six-site unit cell with positions specified as 
ai = (0,0), a 2 = (0,1), a 3 = (V3/2.3/2), a 4 = (V3,l), 
a 5 = (%/3,0), and a 6 = (V3/2, -1/2). The supercell 
lattice vectors are given by Aj/ 2 = (3vo/2, ±3/2). 

Within CDMFT as well as DCA, the interacting lattice 
Green's function in the site basis is defined as 

Gij{iu) n ) = ^ [ iuj n + A 1 - M k ) - E{tU3 n )]T. , (2) 

k 

where u> n = (2n + l)irT are Matsubara frequencies and 
T is the temperature. At half-filling, the chemical po- 
tential is fi = U/2. The k sum extends over the reduced 
Brillouin Zone, /i(k) = — i(k), where t(k) denotes the 
hopping matrix for the superlattice, and £y(iw„) repre- 
sents the self-energy matrix in the site representation. 

Within CDMFT, the elements of t(k) within the unit 
cell given by iy = t for neighboring sites. In addition, 
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hopping between cells yields: 



tl4 
*36 



te 
te 
te 



-ik- Ai 
-ik- A 2 
-ik-A 3 



(3) 



where A3 = A2 — Ai. The hopping matrix t(k) is Her- 
mitian, so that tji — t*j . All other elements vanish. 

To distinguish the hopping matrix elements within 
DCA, we denote them by i#(k). In the real-space ver- 
sion of DCAi£ they are related to those within CDMFT 
via a phase factor: 



U 



U 



— zk-i 



(4) 



where ay = 
lowing matrix elements: 



a,-. This phase relation yields the fol- 
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= h2 
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= te~ ika3i , 





with analogous connections among the Hermitian ele- 
ments. All other matrix elements vanish. 

The cluster Hamiltonian in CDMFT has the hopping 
matrix elements Ek^(k)]y = tfj where tfj = t = 1 for 
first neighbors and ifj = otherwise. In contrast, in DCA 

we find ifj =i= 0.8103 for first and third neighbors and 

H 



tfj — otherwise. 

Within CDMFT as well as DCA 



Gij is a symmetric 



matrix, with site-independent diagonal components Gu. 
In the case of CDMFT, there are three independent off- 
diagonal elements: G12, G13 and G14. Here, Gu, G13 are 
imaginary and G12, G14 are real. Thus, the correspond- 
ing density of states components pn and pn are even 
functions of energy, while P12 and P14 are odd. In the case 
of DCA, translation symmetry is preserved, so that one 
has the additional condition P12 = pu and G12 = G14 
due to the equality of first- and third-neighbor hopping 
interactions t. 

Because of these symmetry properties, it is useful 
to express the lattice Green's function in the diagonal 
molecular-orbital basis whose elements G m (iuj n ) (m = 
1 ... 6) are determined by: 



Gi ; 2 
G3,4 = G$fi 



(Gu 
(Gu 



2G13) i (G14 
G13) ± (G14 - 



H2G 12 ) 
G 12 ). 



(6) 



The unitary transformation Ti m linking the site and 
molecular-orbital bases is defined in Eq. (6) of Ref.@. Ev- 
idently, in CDMFT there are two independent complex 
functions, Gi = — G% and G 3 = — G%. In DCA, the ele- 
ments G m= 3...6 are degenerate and imaginary. The onsite 
and intersite components of the lattice Green's function 
can be derived by inverting Eq. ([6]): 



Gu 
G12 

G13 
G14 



[(Gi 
[(Gi 
[(Gi 
[(Gi 



G 2 ) 
G 2 ) 
G 2 ) 
Gi) 



2(G 3 + G 4 )]/6 
(G 3 - G 4 )]/6 
(G 3 + G 4 )]/6 
2(G 3 -G 4 )]/6. 
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FIG. 1: (Color online) Density of states p(oj) (solid curves) of 
honeycomb lattice and cluster components p m (^) in diagonal 
molecular orbital basis (dashed curves) for (a) CDMFT and 
(b) DCA. For clarity, these components are divided by n c — 6. 
In CDMFT, all density components are non-symmetric and 
orbitals 3 and 4 are doubly degenerate. In DCA, only pi 
and p2 are non-symmetric, while the degenerate components 
p m =3 ...6 are symmetric, uj = defines the Fermi energy for 
half-filling. 



Figure [T] illustrates the uncorrelated density of states 
components in the diagonal molecular orbital basis, 
where p m (w) = -jImG m (w). The total density of 
states is, of course, the same within CDMFT and DCA, 
but its decomposition into molecular-orbital or intersite 
contributions differs for these two schemes. The four 
CDMFT densities shown in panel (a) are non-symmetric 
and satisfy the relations P2(w) = Pi{— u) and pi{uj) = 
Ps(—oj). The corresponding DCA densities are plotted 
in panel (b). In this case, only pi(co) = P2(~ are non- 
symmetric, whereas Pz(uj) — P4(w) are symmetric. 

Figure [5] (a) shows the Brillouin Zone of the honey- 
comb lattice together with the three times smaller re- 
duced Zone. Panel (b) illustrates the contributions to 
the density of states stemming from the outer k regions 
KMK'M' and the inner regions TM'K'. These two con- 
tributions overlap slightly since the point K' does not lie 
half-way between T and M . Thus, the low-energy part of 
the density of states (denoted as K) extends up |<j| < 2, 
while the high-energy part (denoted as T) corresponds to 
the window 1.75 < < 3. The comparison with Fig. Q] 
(b) shows that the diagonal elements of the DCA den- 
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FIG. 2: (Color online) (a) Segment of Brillouin Zone of 
honeycomb lattice (solid red lines). The reduced Zone 
(dashed blue lines) is obtained by folding the Dirac points 
K onto F and the M points onto M' . (b) Decomposition of 
density of states into low-energy contribution (denoted as K) 
corresponding to outer regions KMK'M' and high-energy 
contribution (denoted as T) corresponding to inner regions 
TM'K' of large Brillouin Zone. 



sity of states correspond to the distributions indicated in 
Fig. [2] (b). Thus, pi j2 (w) account for the energy bands 
in the inner regions TM'K' and p^^uj) for those in the 
outer regions KMK'M' of the original Brillouin Zone. 
The momentum regions shown in Fig. [5] (a) therefore 
specify the appropriate tiling of the Brillouin Zone within 
the DCA. 

The self-energy matrices in CDMFT and DCA satisfy 
the same symmetry properties as the lattice Green's func- 
tions so that they can be diagonalized in the same man- 
ner. These diagonal elements will be denoted as £ m (iw n ). 
In the site basis the components En and £13 are imagi- 
nary, whereas £12 and £14 are real. As translation sym- 
metry is not obeyed in CDMFT, £12 and £14 differ, while 
in DCA they coincide. 

We point out that, although the hopping matrix ele- 
ments t(k) in CDMFT and DCA differ only by a unitary 
transformation as indicated in Eq. the same does 
not hold for the respective self-energy matrices. As dis- 
cussed below, the preservation of translation invariance 
in DCA and its absence in CDMFT give rise to funda- 
mentally different physical solutions which severely af- 



fect the phase boundaries. Thus, the DCA and CDMFT 
self-energy matrices are not simply related via a unitary 
transformation. 

Severe differences of this kind between DCA and 
CDMFT do not arise in the case of the Hubbard model 
for the square lattice, where the cluster Hamiltonians 
maintain the same symmetry. The only difference is that 
the hopping interaction between neighbors is changed 
from t = 1 in CDMFT to t = 1.273 in DCA. As a result, 
these cluster schemes lead to a less dramatic reorgani- 
zation of spectral weight among the cluster components 
than in the case of the honeycomb lattice. 



III. EXACT DIAGONALIZATION 

To avoid double-counting of Coulomb interactions in 
the quantum impurity calculation, the self-energy must 
be removed from the six-site cluster in which correla- 
tions are treated explicitly. This removal yields the bath 
Green's function matrix 



G (ioj„) = [G(iu) n ) 1 -I- T,(iu! n )}~ 



(8) 



Within the ED approach, this bath Green's function of 
the infinite lattice is projected onto the corresponding 
function of a supercluster consisting of n c = 6 correlated 
sites within the unit cell plus a bath consisting of rib dis- 
crete levels. Here, we choose n& = 6, so that the total 
number of levels of the supercluster is n s = n c + n& = 12. 
Within the diagonal molecular-orbital basis, this projec- 
tion implies 



G , m (io; n ) « G$ t7n (iu} n ) 

= \iuj n + [x - e T , 



E 

k=l 



<mk I 



,(9) 



where e m =i...6 denotes impurity levels and 6/5=7,, .12 bath 
levels. The bath levels are defined relative to the chemical 
potential. We assume that the molecular orbitals couple 
to independent baths so that the hybridization matrix 
elements are also diagonal in this representation: V m k — 
Sm+6,kVk- Fig. [3^a) illustrates the impurity and bath 
levels in the diagonal molecular orbital basis. Panel (b) 
shows the equivalent representation when the impurity 
orbitals are transformed to the original site basis. The 
bath remains unchanged and the hopping terms in this 
basis are given by Vik — J2 m Ti m V m k- This picture differs 
from the one in which also the bath is treated within the 
site basis (see below). 

To determine the bath levels and hopping terms V m k 
we minimize the difference 

M 

Diff m = W£\G 0>m (iu n ) - G c [ m (ioj n )\ 2 , (10) 

n=0 

where M w 2 10 is the total number of Matsubara points 
and the weight function — l/co^ 1S introduced to 
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FIG. 3: (Color online) (a) Cluster levels in molecular orbital 
basis. There are six independent terms connecting orbital 
levels e m= i...6 (red dots) and bath levels tk=7.. .12 (blue dots) 
via hopping integrals Vfc=7...i2- In CDMFT (for fixed impurity 
levels) one has: £1,2 = ^24, £3,4 = £5,6 = it, and £7 = 

— 68, £9 = £11 = — £10 = — £12, Vt = Vg, Vg — VlO = Vll = 

V12. Thus there are four independent bath parameters. In 
DCA, £1,2 = =F3t, £3. ..6 = £9. ..12 = 0, i.e., there are only 
three independent fit parameters, (b) Cluster levels in site 
basis i = 1...6 (green dots) connected to molecular orbital 
bath levels £fc=7...i2 (blue dots) via hopping integrals Vik- For 
clarity, the hopping interactions between impurity sites are 
not shown. Representations (a) and (b) are equivalent since 
they are connected via the unitary transformation T between 
impurity sites i = 1...6 and orbitals m — 1...6. The bath 
molecular orbital levels in (b) are the same as in (a). Thus, 
although the cluster sites have identical levels at zero energy, 
the bath levels maintain the orbital symmetry. 



give more weight to the low-frequency region. We usually 
take N = 1 or N = 2. Note also that both Green's func- 
tions in the above expression approach l/iui n for large 
oj n . Thus the difference defined in Eq. (jTU)) automati- 
cally focuses on the low-energy region. This is not the 
case when the differences of the inverse Green's functions 
are minimized instead. The reason is that the hybridiza- 
tion functions corresponding to C?o,m and Gq m are not 
normalized to the same asymptotic amplitudes. To start 
the iterative procedure, we use bath parameters obtained 
for the uncorrelated system, or from a converged solution 
for nearby Coulomb energies. The resulting and V m k 
are usually very stable against variations of initial condi- 
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FIG. 4: (Color online) Comparison of lattice bath Green's 
function Go !m (iw„) (solid red curves) and cluster Green's 
function (dashed blue curves) for U — 3 and T = 0.01. As the 
density of states for m = 3 is symmetric in DCA (see Fig. 1 
(b)), Go, 3 is purely imaginary, while Go,i is complex. Thus, 
the latter function is fitted with two parameters, whereas Go, 3 
involves only one fit parameter. The solid and dashed curves 
for Go,i are indistinguishable. 



tions. 

In the CDMFT calculations discussed in Ref. [|1 not 
only the bath levels £& and hopping elements Vk were 
used as parameters in the fit of Go jm (iu> n ), but also the 
impurity levels e m . Since the expression Eq. ([5]) en- 
sures the correct asymptotic behavior, the variation of 
e m yields slightly better accuracy of the fit at the lower 
Matsubara points. For each diagonal component Go. m 
three fit parameters are then available. As there are only 
two independent complex functions Go,m, the total num- 
ber of parameters to fit the bath is six. As shown in 
Fig. 24 of Ref. \M for U = 4 and T = 0.01, this proce- 
dure yields a surprisingly good reproduction of the lattice 
bath Green's function via the cluster Anderson Green's 
function, in spite of the fact that we use only one bath 
level per impurity orbital. The reason for this good fit is 
the semi-metallic nature of the honeycomb lattice, giving 
rise to a vanishing density of states at the Fermi level. 
In contrast, in ordinary correlated metals and the trian- 
gular or square lattice Hubbard models, the density of 
states of the infinite lattice is finite, so that a successful 
fit to a cluster Green's function usually requires at least 
two bath levels per orbital and restriction to not very low 
temperatures (typically T > 0.01). 

In the DCA calculations presented below, we fix the 
impurity levels e m at their nominal cluster values. Thus, 
£1.2 = and £3,4 = 0. The latter value reflects the 
fact that the DCA density of states components p3^(oj) 
are even functions of energy. Thus, the fit of Go,m=i,2 
involves two parameters (the bath level £7 = — e$ and the 
hopping element V7 = Vg), while Go, m= 3,4 includes only 
the hopping element Vg = Vlq as fit parameter. 
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Figure [4] illustrates the quality of the fit of Gq within 
ED DCA for U = 3 and T = 0.01. The parameters used 
in these fits are: a = -3t = -2.4309, e 7 = -1.85694, 
V 7 = 0.26270 for to = 1 and e 3 = e 9 = 0, Vg = 0.86701 for 
m = 3. As pointed out above in the case of CDMFT, the 
excellent representation of the lattice Green's function 
via the cluster Green's function using only one bath level 
per impurity orbital is related to the vanishing density of 
states at the Fermi level. 

The diagonalization of the supercluster Hamiltonian is 
conveniently carried out in the site basis. At low temper- 
atures only few excited states need to be included in the 
evaluation of the cluster Green's function Gfj(ioj n )- The 
diagonalization can then be performed very efficiently by 
making use of the Arnoldi algorithm. Details concern- 
ing this procedure are provided in Refs. @ bmce 
the cluster Green's function obeys the same symmetry 
properties as the lattice Green's function, it is diagonal 
in the molecular-orbital basis. These elements will be 
denoted as G^(iuj n ). The diagonal cluster self-energy 
components are then given by an expression analogous 
to Eq. ©: 



(11) 



The key physical assumption in DMFT is now that this 
cluster self-energy provides an accurate representation of 
the lattice self-energy. Thus, 



S m (iw„) w S^(iw„). 



(12) 



In the next iteration, these self-energy components are 
used as input in the lattice Green's function Eq. ©. 
In the diagonal molecular-orbital basis the DCA lattice 
Green's function is given by 

G m (iu n ) = ^ [ ioj n + M - f^h^f - T,(ioj n )] mm . 

(13) 

We note here that, at real energies, the cluster quantities 
Gm, G§ m and have discrete spectra, while the cor- 
responding lattice spectra associated with the quantities 
G m , Go, m and E m are continuous. 

We close this section by pointing out that we believe 
the projection of the bath Green's function within the 
diagonal molecular-orbital basis discussed above to be 
more general and more flexible than analogous projec- 
tions within the nondiagonal site basis. As mentioned 
above, within CDMFT there are two independent com- 
plex functions Go,m (with nonsymmetric spectral distri- 
butions) that are fitted each with one bath level tk and 
one hopping term Vk (assuming the impurity level e m to 
be fixed). Thus, there are altogether four fit parameters. 
This should be compared to only one fit parameter if 
the site basis is used instead. For symmetry reasons all 
bath levels then are zero so that only the site indepen- 
dent impurity bath hopping element remains as a single 
fit parameter. Introducing a hopping interaction among 
bath levels as was done in Ref. Illl increases the number 
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FIG. 5: (Color online) Comparison of excitation gaps as func- 
tions of Coulomb interaction derived using several cluster 
methods and impurity solvers for paramagnetic phase of hon- 
eycomb lattice: Meng et al: large-scale QMC— , Wu et al.: 
CTQMC CDMFT—, Liebsch: ED CDMFT—, He et al: ED 
CDMFT—, Seki et al: ED VCA 14 . In contrast, both ED and 
CTQMC DCA yield semi-metallic behavior with A = for 
U < 6 (see text). 



of fit parameters from one to two. Actually, since the 
bath can always be represented in a diagonal form, hop- 
ping among bath levels is implicitly included in the diag- 
onal molecular orbital picture with four fit parameters. 
Analogous considerations hold for DCA. Nevertheless, as 
will be shown in the next section, these slightly different 
implementations of ED all yield consistent answers con- 
cerning the variation of the excitation gap as a function 
of Coulomb energy. 



IV. RESULTS AND DISCUSSION 

Figure [5] shows the comparison of the excitation 
gaps obtained for various cluster methods and impurity 
solvers. Near U » 5, all calculations (except DCA, see 
below) predict a Mott phase with a gap A 0.5 — 0.9. 
At U < 4, the CDMFT and VCA results that do not pre- 
serve translation symmetry exhibit a gap tail that per- 
sists down to U — > 0. The differences between these 
results are partly caused by the different temperatures 
used in these studies. In particular, the gap closing near 
U = 3.8 obtained within CDMFT by Wu et alA seems 
to be related to the rather high temperature, T = 0.05, 
employed in the CTQMC calculation. Since the CTQMC 
self-energy agrees well with the ED results, CTQMC 
CDMFT presumably would also yield a gap at lower T. 
Also, the ED calculations in Ref. |6 were carried out at 
T = 0.005, while those in Refs. [Ill HI essentially corre- 
spond to the T — > limit. 

In striking contrast to CDMFT, the translation invari- 
ance of DCA ensures the existence of a semi-metallic 
phase at low values of U . On the other hand, the con- 
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FIG. 6: (Color online) Density of states Aii(u>) = 
— ^lmGn(uj) of honeycomb lattice for several Coulomb en- 
ergies at. Red solid curves: U — 6, dashed curves: U = 3 — 5. 
(a) ED DCA (T = 0.01) (b) CTQMC DCA (T = 0.025). 
For illustrative purpose, only the low-energy range of the ED 
spectra is shown. The dotted curve denotes the bare density 
of states. 



dition S12 = £14 cannot generally be correct for the 
short-range correlations within the unit cell. Thus, at 
Coulomb energies, where local Mott physics dominates 
and long-range translational invariance becomes less im- 
portant, DCA should be less appropriate than CDMFT. 
Indeed, both ED and CTQMC DCA results suggest that 
the semi-metallic phase with A = extends to U > 6, 
i.e., beyond the critical Coulomb energy U c « 3.9 — 4.3 
of the antiferromagnetic phase*^ 

This is illustrated in Fig.[6j which shows the interacting 
density of states obtained in ED and CTQMC DCA for 
several Coulomb energies. The ED spectra were obtained 
by making use of the extrapolation routine ratint^ while 
the CTQMC spectra were derived via the maximum en- 
tropy methodi^i For details concerning the CTQMC cal- 
culations, see Ref. 0. The main effect of Coulomb in- 
teractions is seen to be the usual band narrowing and 



FIG. 7: (Color online) Green's function components G\i(ioj n ) 
(i — 1,2,3) of honeycomb lattice as functions of Matsubara 
frequency calculated within (a) ED DCA and (b) CTQMC 
DCA at T = 0.01. Solid red curves: (7 = 4; dashed curves: 
(7=1-3. 



effective mass enhancement, as found in weakly corre- 
lated systems. In contrast, the corresponding ED and 
CTQMC CDMFT spectra for U = 5 reveal a large Mott 
gap of about A = 0.6 (see Fig. \5$M 

The persistence of semi-metallic behavior at large U 
within DCA is related to the fact that the enforcement of 
translation symmetry is achieved at the expense of equat- 
ing first- and third-neighbor interactions in the cluster 
Hamiltonian. The self-energy in the site basis then sat- 
isfies the condition £12 = £14, whereas in CDMFT £14 
is noticeably smaller than £ 12 iSiifl 

The good correspondence between the DCA spectra 
obtained within ED and CTQMC is a consequence of 
the nearly quantitative agreement of the lattice Green's 
functions Gu(iuj n ) which are shown in Fig. [7] As pointed 
out in the preceding section, for symmetry reasons Gn 
and G 13 are imaginary, while G12 = G14 are real. 
Both impurity solvers yield ImGii(m> n ) — > in the 
limit w„ — ¥ 0, implying that the local density of states, 
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FIG. 8: (Color online) Green's function components Gu(iu)„) 
(i — 1,2,4) of honeycomb lattice as functions of Matsubara 
frequency calculated within ED CDMFT. Solid red curves: 
U — 4: dashed curves: (7 = 1 — 3. 



p(ui) = — ilmGn(w) vanishes at oj = 0. Also, both 
schemes indicate that with increasing values of U the 
initial slope of ImGn and ImGi3 increases. Thus, the 
Dirac cones become steeper and spectral weight is shifted 
towards the Fermi level. 

The results obtained within DCA differ in two qualita- 
tive aspects from those derived previously in CDMFT. 
As shown in Fig. [31 the components G12 and G14 in 
CDMFT do not coincide. Moreover, the initial slopes 
of ImGn and ImGi3 become smaller with increasing U 
rather than larger as within DCA. In Ref.^it was demon- 
strated that for U > 4 a Mott gap opens in the density 
of states, in approximate agreement with the large-scale 
QMC calculations by Meng et alX At smaller values of 
U, a tiny gap or pseudogap was also found (see below), 
which is however difficult to resolve within ED at finite T. 
As the opening of a gap in the density of states implies 
a reduction of |Im Gn(iw n )| at small values of w„, the 
results shown in Figs. [7] and [5] underline the fundamental 
difference between DCA and CDMFT for the honeycomb 
lattice: Whereas DCA yields a weakly correlated semi- 
metal, CDMFT gives rise to insulating behavior even at 
small U . 

To illustrate the effect of Coulomb correlations in more 
detail, we show in Fig. [9] the self-energy components in 
the site basis for several values of U . The correspond- 
ing results obtained within CTQMC DCA are depicted 
in Fig. 1101 There is good overall correspondence between 
these two impurity solvers, except for slightly different 
magnitudes of the off-diagonal components. We note, 
however, that Re £12 and ImEi3 are approximately one 
and two orders of magnitude smaller than ImEn, re- 
spectively. As can be seen in Fig. [JJ these differences 
have only a minor effect on the variation of the Green's 
function components with increasing Coulomb energy. 

The crucial question in the case of the honeycomb lat- 
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FIG. 9: (Color online) Self-energy components Ti\i{iu n ) (i — 
1,2,3) of honeycomb lattice as functions of Matsubara fre- 
quency calculated within ED DCA for U = 1 - 4 at T = 0.01. 



tice is how Coulomb correlations influence the energy 
bands in the vicinity of the Dirac points. The self-energy 
at these points can be shown to have the simple form: 6 



T,(K, iuj n ) 1=3 itu n a 



iu n (l - a) 
where the coefficients are given by 



u) n 0, 



a = Im [Su(iw„) - i: 13 (iuj n )]/u! ri 
b = Re [Si 2 (iw„) - Ei 4 (iw n )] 



(14) 



(15) 
(16) 



in the limit uj n — > 0. Thus, T,(K,iuj n ) is imaginary 
as expected for particle-hole symmetry near the Dirac 
points. Moreover, this self-energy consists of metal- 
lic (~ iw n ) and insulating (~ l/itu n ) contributions, 
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FIG. 10: (Color online) Sell-energy components Eii(iw„) 
(i = 1,2,3) of honeycomb lattice as functions of Matsubara 
frequency calculated within CTQMC DCA for U = 1 - 4 at 
T = 0.01. 



where the latter term is a direct consequence of the 
fact that S12 ^ £14. The presence of this term im- 
plies ReT<(K,uj) ss b 2 /[uj(l — a)} at real w. In the low- 
temperature limit, this expression yields an excitation 
gap of magnitude A 1=3 2-^/jcf, where c = b 2 /(l — a). A 
similar insulating contribution to the self-energy was re- 
cently found in Ref. Presumably, this insulating term 
is also present in the ED calculations reported in Refs. El 
and iTJ. At finite T, the gap is smoothened out so that 
it becomes difficult to determine its boundaries. In con- 
trast, as discussed in Section II, DCA preserves the bulk 
symmetry, so that £12 = £14 and A = 0. Thus, the DCA 
self-energy at the Dirac points is purely metallic, where 
the increasing magnitude of the coefficient a implies in- 



creasing quasi-particle broadening and shift of spectral 
weight towards the Fermi level as U increases. From the 
initial slope of Im£n at U = 4 we obtain an effective 
mass enhancement of about m*/m ~ 1.3. 

The above discussion demonstrates that the presence 
or absence of the insulating contribution to Y,{K, iuj n ) is 
not caused by the impurity solver used in the CDMFT or 
DCA calculations. In fact, the good agreement between 
ED and CTQMC, for both CDMFT and DCA, suggests 
that in the case of the honeycomb lattice one bath level 
per impurity orbital is sufficient for an accurate fit of the 
bath Green's function. The reason is that, because of 
the semi-metallic nature of the honeycomb lattice, the 
projection of the bath Green's function of the infinite 
lattice onto a finite-cluster Anderson Green's function is 
not plagued by the low-energy-low-temperature discrep- 
ancies that usually occur in the case of correlated metals. 
In these systems at least two bath levels per impurity or- 
bital are typically required and very low temperatures 
must be avoided*^ 



V. SUMMARY 

The role of Coulomb correlations in Hubbard model 
for the honeycomb lattice has been studied within finite- 
temperature exact diagonalization and continuous-time 
quantum Monte Carlo combined with the dynamical clus- 
ter approximation. The unique feature of DCA is that it 
preserves the translation invariance so that the system at 
small values of U is semi-metallic. In contrast, CDMFT 
violates translation symmetry which implies the opening 
of an excitation gap at arbitrarily small U , regardless 
of the impurity solver. This gap is therefore an artifact 
caused by the lack of long-range crystal symmetry and 
does not correspond to a true Mott gap. At larger values 
of U, however, many-body interactions are dominated by 
short-range correlations and translation symmetry seizes 
to be important. DCA then becomes less accurate since it 
overemphasizes semi-metallic behavior. Thus, for U w 5, 
CDMFT is preferable and reveals a Mott gap in qualita- 
tive agreement with large-scale QMC calculations. 

In the case of the honeycomb lattice, DCA and 
CDMFT may therefore be viewed as complementary clus- 
ter approaches. As DCA preserves translation symme- 
try, it is more appropriate in the semi-metallic phase at 
small U where long-range order is a prerequisite for the 
description of the weakly correlated Dirac cones. The 
condition £12 = £14 which guaranties this symmetry, 
however, is unrealistic at larger U, when short-range cor- 
relations within the six-site unit cell begin to dominate. 
Thus, in the region of the Mott phase, CDMFT is more 
suitable. As a result of these inherent limitations of both 
cluster schemes, the critical Coulomb interaction defining 
the precise boundary between these phases is at present 
difficult to determine within either CDMFT or DCA. We 
emphasize that this difficulty is not related to the finite 
size or symmetry of the bath used in ED. On the contrary, 
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within CDMFT as well as DCA, the ED self-energies 
agree well with the corresponding CTQMC results. 

It is interesting to inquire why the remarkable differ- 
ence between CDMFT and DCA for the honeycomb lat- 
tice discussed in this paper does not also manifest itself 
in other systems, such as the Hubbard models for square 
and triangular lattices. In these cases, long-range order is 
mainly responsible for the logarithmic divergence of the 
van Hove singularities of the density of states. Thus, any 
lack of perfect translation symmetry would give rise to a 
rounding of this peak, an effect that would be difficult to 
distinguish from broadening induced by finite tempera- 
ture and quasi-particle damping. In contrast, any round- 



ing of Dirac cones induces the opening of a gap. In this 
regard, the Dirac cones of the honeycomb lattice corre- 
spond to a rather peculiar special situation that does not 
arise in most cases which have been studied previously 
within CDMFT and DCA at finite temperatures. 
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